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ABSTRACT: 

The glycolytic enzyme pyruvate kinase (PK) is required for cancer development, 
and has been implicated in the metabolic transition from oxidative to fermentative 
metabolism, the Warburg effect. However, the global metabolic response that follows 
changes in PK activity is not yet fully understood. Using shotgun proteomics, we 
identified 31 yeast proteins that were regulated in a PK-dependent manner. Selective 
reaction monitoring confirmed that their expression was dependent on PK isoform, 
level and activity. Most of the PK targets were amino acid metabolizing enzymes 
or factors of protein translation, indicating that PK plays a global regulatory role in 
biosynthethic amino acid metabolism. Indeed, we found strongly altered amino acid 
profiles when PK levels were changed. Low PK levels increased the cellular glutamine 
and glutamate concentrations, but decreased the levels of seven amino acids including 
serine and histidine. To test for evolutionary conservation of this PK function, we 
quantified orthologues of the identified PK targets in thyroid follicular adenoma, a 
tumor characterized by high PK levels and low respiratory activity. Aminopeptidase 
AAP-1 and serine hydroxymethyltransferase SHMT1 both showed PKM2- concentration 
dependence, and were upregulated in the tumor. Thus, PK expression levels and 
activity were important for maintaining cellular amino acid homeostasis. Mediating 
between energy production, ROS clearance and amino acid biosynthesis, PK thus plays 
a central regulatory role in the metabolism of proliferating cells. 



INTRODUCTION 

Cancer cell growth is characterized by re- 
configuration of cellular metabolism. Different routes 
of intermediary and biosynthetic metabolism, including 
nucleic acid synthesis, amino acid biosynthesis, and 
energy production are reconfigured [1-4]. A characteristic 



feature of many tumors is an augmentation of non- 
oxidative (fermentative) metabolism at the expense of 
oxidative metabolism (respiration), although the latter 
would be more efficient for energy production. Known 
as the 'Warburg effect', this metabolic reconfiguration is 
considered as a 'hallmark' of cancer, and may have high 
therapeutic potential [5, 6]. 
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Recent results indicate that the Warburg effect could 
be the result of metabolic balancing processes which avoid 
the accumulation of toxic reactive oxygen species (ROS) 
maintaining metabolic homeostasis [7, 8]. Indeed, cancer 
cells suffer from high loads of ROS which may make 
them vulnerable in pro-oxidant therapies [8-10]. ROS 
originate from biochemical reactions, i.e. the degradation 
of fatty acids, and leakage from complex I and III of 
the respiratory chain [11, 12]. Clearing of these reactive 
molecules prevents oxidative stress, and is conducted 
by large machinery which involves central carbon 
metabolism as essential metabolic component. Fluxes of 
glycolysis and the adjacent pentose phosphate pathway 
(PPP) are adjusted when cells encounter an oxidative 
environment [13]. This metabolic reconfiguration fuels 
the antioxidative machinery with its main redox cofactor 
NADPH, required by the glutathione, peroxiredoxin and 
thioredoxin systems [14, 15], and is involved in regulating 
gene expression during oxidative conditions [16]. 

The metabolic enzyme pyruvate kinase (PK) is 
required for coordination of these pathways when cells 
encounter increased ROS levels. Lowering PK activity 
in yeast and human increases resistance to oxidants [17, 
18]. At least in part, this adaptation is mediated through 
accumulation of the PK substrate phosphoenolpyruvate. 
PEP causes feedback inhibition of several glycolytic 
enzymes, including triosephosphate isomerase (TPI) 
[17]. Lowered TPI activity increases oxidant resistance 
as studied in yeast and C. elegans, and acts through re- 
directing the metabolic flux from glycolysis into the PPP. 
This increases the NADP + reduction rate [13], and triggers 
an adaptation of gene expression towards oxidative 
conditions [16]. 

In this study, we focused on the regulatory events 
upstream and downstream of PK, picturing PK induced 
proteome changes in yeast and human. For this, we 
combined workflows of quantitative proteomics which 
circumvent the use of isotope labeled standards ('label- 
free proteomics'). In yeast, low levels of PK were 
sufficient to induce changes in the biosynthetic proteome, 
primarily affecting protein biosynthesis and amino acid 
metabolism. Subsequently, ion exchange chromatography 
was used to profile free amino acid levels in yeast strains 
possessing different PK activities. Low PK activity caused 
a decline in the concentration of arginine, aspartic acid, 
histidine, lysine, serine, threonine, and valine, and an 
increase in glutamine and glutamate, thus triggered a 
broad reconfiguration of the yeast amino acid profile. 

In human thyroid follicular adenoma, we observed 
upregulation of two orthologues of the identified 
yeast proteins, aminopeptidase AAP-1 and serine 
hydroxymethyltransferase SHMT 1 , and provide evidence 
that their expression level might quantitatively correlate 
with PKM2 expression in the cancer state. Thus, in yeast 
and humans, PK emerges as regulator of carbohydrate, 
energy and also amino acid metabolism. 



RESULTS 



PK levels target cellular biosynthetic metabolism 

Although of reduced catalytic activity [18, 19], 
the expression level of the predominant human pyruvate 
kinase isoform PKM2 is increased in cancers [20]. In 
yeast, changes in PK isoform and expression levels 
are observed when cells switch from fermentation to 
oxidative metabolism [17, 21], indicting that PK plays an 
evolutionary conserved role in the underlining metabolic 
reconfiguration. 

Previously, we have shown that lowering PK 
activity in yeast increases the rate of respiration, and at the 
same time augments oxidant resistance [17]. This indicates 
this enzyme is required for coordination of respiration 
and ROS clearance. Here, we used mass spectrometry 
based proteomics to picture the global physiological 
consequences triggered by a change in PK levels. We 
compared two engineered yeast strains that just differ 
in the promoter (TEF1 for high level expression and 
CYC1 for low level expression) used to express yeast 
Pyklp",' the high active yeast PK isoform [17, 21]. MS/ 
MS shotgun proteome profiles of these two yeast strains 
were recorded using a nanoLC (Ultimate 3000 (Dionex)) 
coupled QqTOF mass spectrometer TripleTOF5600 
(AB/Sciex) [22]. Obtained MS/MS spectra were used 
to identify peptides and proteins by a database search 
with ProteinPilot™ software (AB SCIEX). A typical run 
yielded the identification of 922 proteins at global FDR of 
<1% from the UniprotKB/Swiss-Prot database 57.15. At a 
global FDR of 1% on the peptide level 6393 peptides and 
on spectrum level 72256 spectra were identified. 

The identified proteins and peptides were analyzed 
with PeakView™ software (AB SCIEX), where mass 
and extracted ion chromatograms (XICs) were obtained 
for the confidentially identified peptides (Peptide ID 
with unique peptides at a global false discovery rate < 
1%). The same peptides were extracted for all technical 
and biological replicates. The extraction window for all 
peptide XICs was set to 0.01 Da to increase specificity 
and improve the quantification results. Typical elution 
profiles for peptides show a width of 20 seconds (FWHM) 
at the HPLC conditions selected to ensure a number of > 
8 points for quantification, typically at <35 points across 
peak base. Extracted ion chromatograms of the tryptic 
Pyklp peptide K.TNNPETLVAL.R of cells with different 
PYK1 expression levels, 19 data points per peak, illustrate 
clear separation between high (A) and low (B) level 
expression of PK in three replicates each (Fig 1A). For 
relative quantification, sum of the total peak area was used 
to normalize peak areas between the different samples. Fig 
IB illustrates the normalized peak area for all identified 
Pyklp-belonging peptides over technical and biological 
replicates. In total, 42.6% of identified proteins met the 
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Figure 1: Label-free shotgun proteome profiles of yeast expressing Pyklp at high and low level. (A)Tryptic digests 
generated from yeast protein extracts were analyzed on an AB SCIEX QqTOF mass spectrometer. Extracted ion chromatograms for the 
Pyklp tryptic peptide K.TNNPETLVAL.R from three LC-MS/MS shotgun runs each of yeast expressing Pyklp at high level (,A' samples) 
or low level (,B' samples) .(B)Peak area of all Pyklp-associated peptides in cells expressing Pyklp at and high and low level, normalized 
to the total peak area of the respective run. Replicate injections of the same sample (LC-MS technical replicates) are marked with the same 
label; Numbers (1) & (2), respectively (3) & (4) indicate sample preparations and tryptic digests from independently grown yeast cultures. 
(C)Principal component analysis after pareto scaling of the shotgun proteome profiles, generated with Marker View (AB/Sciex). A change 
in the Pyklp expression level clearly separates the proteomes. Determined by two tailed T-test, 31 proteins (Table 1) were significantly (p 
< 0.05) regulated and responsible for the separation of the proteome profiles. 
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stringent criteria of being identified with > 3 unique PK targets respond to level changes of both yeast 

peptides in all 19 data-dependent acquisition runs, thus PYK isozymes 

were used for quantification. 

We first analyzed the shotgun proteome profiles by 
principal component analysis (PCA) after pareto sorting. 
Profiles obtained from both from technical and biological 
replicates showed clear clustering according to whether 
Pykl was expressed at high or low level (Fig 1C). Thus, 
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a change in the expression level of the Pyklp enzyme 
caused a reproducible change in the yeast proteome. 

We continued our investigations by an in depth 
verification of the proteomic results making use of a 
quantitative MS technique, selective reaction monitoring 
(SRM). In difference to shotgun proteomics, SRM assays 
are generated in a targeted way and on the base of the 
protein/peptide sequence or its known fragmentation 
spectrum [23]. With SRM, a much lower number of 
molecules per sample can be quantified (here we choose 
nine proteins randomly), but these assays are more 
sensitive due to improved signal to noise ratio, therefore 
facilitating faster runtimes, and are more reproducible, 
and are thus used in diagnostics [23-25]. We established 
and validated SRM assays for a representative number of 
proteins which were significantly (p < 0.05) responsible 
for separating the proteomes in the PCA (Fig 1 C, Table 



1). The SRM assays were applied to BY4741, the two 
Pyklp strains used in the QqTOF analysis (Fig 1A), and 
two further yeast strains (TEF -Pyk2 and CYC -Pyk2) 
which express the low active Pyk2 isozyme at high and 
low level (please see Fig 2A for an overview, strains were 
in detail described in [17]) 

Label-free SRM confirmed the PYK-dependent 
regulation as identified with the shotgun platform for 
7 of the 9 proteins (Fig 2, Suppl. Table 1 for numeric 
values). Interestingly, all Pyklp regulated proteins were 
also dependent on the Pyk2p expression level, thus 
they were not regulated in an isoform specific manner. 
Remarkably, we observed that the expression pattern 
grouped the PK targets according to their biological 
function. The first group of co-regulated proteins was 
formed by the glycolytic enzymes glycerol-3 -phosphatase 
dehydrogenase, phosphoglycerate mutase, and enolase 2. 
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Figure 2: Targeted quantification of PK regulated genes by selected reaction monitoring. (A)Overview of the yeast model 
strains [17] that differ in overall PYK activity, Pyk isoform, and expression level. (Left bars) Total PK activity changes over BY4741 > 
TEFpr-PYKl > TEFpr-PYK2 > CYCpr-PYKl to CYCpr-PYK2, as described previously [17] . (Right bars). High expression of Pyklp 
and Pyk2p under control of the TEF1 promoter, low expression under control of the CYC1 promoter [17]. (B)Glycolytic enzymes 
glycerol-3-phosphatase, phosphoglycerate mutase 1, and enolase 2 are upregulated in yeast strains with lower PK activity. (C)Ribosomal 
proteins quantified using single reaction monitoring (60S acidic ribosomal protein P2-alpha and 40S ribosomal protein S9-A and S9-B) 
are downregulated in yeast strains with lower PK activity. (D)PYK targets Sari and Hsp82 are regulated in strains with altered PYK1 
expression, but are not correlating with the overall enzyme activity. (A-D) Determined expression values are given in the Suppl. Table 1 . 



www.impactjournals.com/oncotarget 



1360 



Oncotarget2012; 3: 1356-1369 




Amino acid & secondary metabolism 
Translation & ribosome biogenesis 
Chaperones 

Glycolysis & pentose phosphate 
pathway 

Nucleotide biosynthesis 
Other 



Gene 


Function 


Acc. No 


P" 

value 


Human 
orthologue 


Acc. No 




AAP1 


Arginine/alanine aminopeptidase 


P37898 


3.00E-008 


AAP-S (PSA1) 


P55786 




SHM1 


Serine hydroxymethyltransferase 


P37292 


2.83E-005 


SHMT1 


NP_004160 




PSA1 


GDP-mannose pyrophosphorylase 


P41940 


5.55E-005 


GMPPB 


NP 037466 

2 




ASN1 


Asparagine synthetase 


P49089 


4.12E-005 


ASNS 


P08243 




SAM2 


S-adenosylmethionine synthetase 


P19358 


3.40E-005 


METK1 


Q00266 




HIS4 


Histidine bio synthesis 


P00815 


2.17E-011 








MET6 


Methionine synthase 


P05694 


8.71E-007 









B) amino acids, decrease in concentration 

histidine threonine 



lysine 



i-l 120% l 120% l 

Ilk llii Ilk 

^yS* <r <T <T jsa* <r <T <T oss* <T <r <T 



4* 4' (jy- r & 

C) amino acids, increase in concentration D) 

glutamine glutamate 

200%i ** 160%-, ,» 

kl kill 



.rf- / ^ ^ 

# 4 i 4* £ 



r 4 4 4' <f* 



/ 4* 4' <?' 



Significantly altered amino acids 






(p < 0.05) 




decreased 






Arginine 


Aspartic acid 




Histidine 


Lysine 




Serine 


Threonine 




Valine 


increased 


Glutamine 


Glutamic acid 





Figure 3: Amino acid and biosynthethic metabolism are targets of yeast PK. (A)PK targets identified in by proteome profiling 
belong to primary metabolism as well to amino acid- and protein biosynthesis (upper panel). The table lists the genes which play a role 
in amino acid and protein biosynthethic metabolism, and their direct human orthologue (if applicable) (B)Pyruvate kinase affects the 
concentration of nine amino acids. Illustrated is the concentration of histidine, threonine and lysine (decreased) or (C)Glutamine, glutamate 
(increased) relative to their level in BY4741 (wild-type). Significant (p<0.05, two-tailed T test) concentration changes are marked (**). (D) 
In total, nine out of 17 amino acids were found to be significantly regulated (p< 0.05, two-tailed T test), by PK activity; these amino acids 
except glutamine and glutamate were present in lower concentrations in yeast with low PK activity. 
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All these proteins were gradually upregulated as lower 
overall activity of PK was (Fig 2B). Second, ribosomal 
proteins grouped together, and showed the exact opposite 
trend being down-regulated in strains with lower PK 
activity (Fig 2C). Finally, there was the third group of 
proteins which was clearly affected by altered PYK 
expression, but did not correlate with the enzymatic 
activity (Fig 2D). This group was represented by SARI, 
a small COPII coat GTPase and HSC82/ HSP82, an ATP- 
dependent molecular chaperone. These proteins were 
expressed at quite similar levels in TEF-PYK1 and TEF- 
PYK2 yeast, at higher levels in CYC-PYK1 and CYC- 
PYK2 yeast, and even higher levels in the wild-type strain 
BY4741(Fig 2D). 

Hence, low PK activity led to an up-regulation of 
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glycolytic proteins assayed, and at the same time to a 
downregulation of tested components of the translational 
machinery. Furthermore, the proteome study identified 
PK target proteins which responded in an activity 
independent manner. Finally, all proteins found to be 
regulated by PYK1 were also controlled by PYK2, Thus, 
both expression level and enzymatic activity of PK show 
to have an effect on protein expression regulation. 

PK expression levels affect the concentration of 
free amino acids 

Virtually all proteins responsible for the clear 
division of the PYK proteomes were involved in energy, 
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Figure 4: AAP-S and SHMT1 are upregulated in thyroid follicular adenoma and correlate with the PKM2 expression 
level. (A)PKM2, but not PKM1 concentration, distinguishes thyroid follicular adenoma from normal thyroid. PKM1 and PKM2 were 
quantified in tumor and the controls [20], the absolute PKM concentrations in individual biopsies are plotted. (B)Activity of the respiratory 
chain is reduced in thyroid follicular adenoma. Enzyme activity of complexes I - V of the four tumors and control tissues relative to the 
activity of citrate synthase. (C)Puromycin sensitive aminopeptidase (AAP-S) is upregulated in thyroid follicular adenoma, and correlates 
with PKM2 in the tumor state. Relative level of AAP-S in tumor and normal thyroid (left panel), AAP-S expression level plotted versus the 
PKM2 concentration in the same biopsy (right panel). (D)Serine hydroxymethyltransferase 1 (SHMT1) shows clear trends towards higher 
expression in thyroid follicular adenoma, and seem to correlate with PKM2 in the tumor state. Relative level of SHMTlin tumor and 
normal thyroid (left panel), SHMTlexpression level plotted versus the PKM2 concentration (right panel). (A-D) p- values were determined 
with two-tailed T test, n.s:. non significant 
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biosynthetic and intermediary metabolism as well as 
translation (including ribosome biogenesis and protein 
folding) (Fig 3A). This indicates that changing PK levels 
triggered metabolic reactions and pathways required for 
the biosynthesis of proteins (biosynthethic metabolism). 
Importantly, the list of significantly regulated yeast genes 
contained several enzymes involved in the synthesis of 
amino acids. Amino acid concentrations are central for 
biosynthetic metabolism, important for growth control, 
and at the same time implicated in cancer metabolism [1, 
26, 27]. Therefore, we continued our survey by quantifying 
17 amino acid levels using ninhydrin detection/ion 
exchange chromatography in wild-type as well as the 
yeast cells with altered Pyklp and Pyk2 expression levels. 
Indeed, PK expression had large effects on the amino 
acid profiles, nine of the 17 profiled amino acids were 
significantly (P < 0.05) altered by a change in the PK 
expression level (Fig 3). Histidine, threonine or lysine, for 
instance were strongly decreased in cells with low PYK 
activity (Fig 3B), whereas glutamine and glutamate were 
detected at higher concentration (Fig 3C). Only eight 
amino acids, including methionine, glycine or valine did 
not significantly respond to a change in the activity of PK 
(Suppl. Table 2). Thus, PK levels influence the expression 
of proteins involved in biosynthetic metabolism, and cause 
a strong signature in the amino acid profile. 

PKM2- dependent regulation of AAP-S and 
SHMP in Thyroid Follicular Adenoma 

Both, the enzymatic/ glycolytic functions of PK, and 
their role in the regulation of redox processes are highly 
conserved from yeast to humans [7, 8]. Therefore we 
speculated that also this newly discovered aspect of the PK 
function could be conserved. Therefore, we selected PK 
targets for which we could identify a human orthologue, 
and tested if they were expressed at sufficient level to be 
quantified by SRM in thyroid follicular adenoma. 

Thyroid follicular adenomas are benign human 
cancers which exhibit clear up-regulation of PKM2 
when compared to normal thyroid tissue [20] (Fig 4A 
for a summary). First, we determined the activity of 
the five respiratory chain complexes in these samples. 
The enzymatic activities were normalized to the citrate 
synthase, which is a marker of mitochondrial mass. A clear 
reduction in respiratory chain activity (complexes I - IV) 
was observed in thyroid follicular adenomas compared to 
normal adjacent thyroid tissue (Fig 4 B). Thus, exhibiting 
upregulation of PKM2, and low activity of the respiratory 
chain, follicular thyroid adenomas represented an 
archetypal tissue exhibiting the Warburg effect. 

Using the QTRAP mass spectrometer, we 
tested if human orthologues of the yeast proteins can 
be quantified by label free LC-SRM in the trypsin 
digested cell extracts [25]. We could identify peptides 



and SRM transitions suitable for robust relative 
quantification of the aminopeptidase AAP-S and serine 
hydroxymethyltransferase SHMT1 (Fig 4). Specific SRM 
signals corresponding to the other three proteins were of 
too low signal or specificity for a reliable quantification 
with our method (data not shown). In adenoma, both 
AAP-S and SHMT1 were present in higher concentration 
compared to normal thyroid tissue (Fig 4C), showing 
that these proteins are upregulated in the tumor state. 
The effect was prominent for AAP-S; the AAP-S levels 
clearly distinguished normal and adenoma tissue (p < 
0.01). Next, we plotted the relative expression level of 
AAP-S versus the concentration of PKM2 (Fig 4C, right 
panel) Both, AAP-S and SHMT1 did not correlate with the 
expression level of PKM2 in the healthy tissue. However, 
both proteins showed correlation with the PKM2 level in 
the tumor state (Fig AC). AAP-S levels exhibited an almost 
perfect linear correlation with the PKM2 concentration. 
Thus, direct orthologues of two yeast PK target genes 
were upregulated in human thyroid follicular adenomas. 
Moreover, their expression correlates with the PKM2 level 
in the tumor, but not in healthy tissue. 

DISCUSSION 

Reaction cascades of central metabolic pathways are 
well characterized, yet still relatively little is known about 
their regulation and coordination. However, homeostasis 
of cellular metabolites and metabolic fluxes is essential 
for cellular growth, and needs to be adapted to the growth 
conditions of all living cells [3, 28-32] Recently, it has 
become clear that the metabolic enzyme PK is required 
for proper coordination of energy- and redox metabolism. 
In its glycolytic role, PK converts phosphoenolpyruvate 
(PEP) to pyruvate, yielding a molecule of ATP, and thus, 
this reaction is required for the net ATP production in 
glycolysis [33]. Interest in PKM2 had increased markedly 
after it had been suggested that PKM2 expression in adults 
could be specific to cancer cells, and that a PKM1 to 
PKM2 isoform switch might be responsible for increased 
fermentative metabolism (Warburg effect) of cancer cells 
[34, 35]. Cancer specificity of PKM2 could however 
not be confirmed in recent studies, as PKM2 was found 
to dominate over PKM1 in most adult tissue, including 
kidney, lung, liver, thyroid, bladder and colon [20], except 
however muscle [34]. Despite not specific to cancer, 
several studies established an essential role of PKM2 in 
metabolic regulation and the oxidative stress response, 
processes of high importance for cancer development 
[17-19, 36-40]. First, the literature appears consistent 
that PKM2 is upregulated on transcript or protein level 
in several malignant and benign tumors. These studies 
included human melanoma [41] childhood leukemia [42], 
renal cell carcinoma, bladder carcinoma, hepatocellular 
carcinoma, colorectal carcinoma, lung carcinoma, thyroid 
follicular adenoma [20], and mouse kidney tumors [36]. 
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Several cancer tissues also expressed the PKM1 isoform, 
however at much lower levels than PKM2. In some but 
not all tumors, also the PKM1 isoform was found to be 
upregulated [20]. 

Conversely, at the same time it has been reported 
that the PK activity is reduced by posttranslational 
modifications in cancer cells, and that this reduction is 
required for their proliferation [18, 19, 43]. Reduced PK 
activity was attributed to PKM2 phosphorylation [19], 
thiol oxidation [18], and acetylation [43]. This indicates 
that PKM2 is upregulated in cancer to increase the protein 
content, but not to increase the overall PK enzymatic 
activity. However, it has also to be noted that it is yet 
unclear to which extent the identified modifications might 
also be found on PKM2 in normal human tissue, and thus, 
and as to which extend they are cancer cell specific. 

Here, we followed the proteome wide consequences 
of the Pyklp downregulation by creating shotgun 
proteome of yeast expressing Pyklp at high and low levels 
using a latest-generation QqTOF proteomics platform, 
a TripleTOF™ 5600 (AB SCIEX) mass spectrometer 
coupled to a high-pressure nanoLC system was used for a 
label free proteome profiling. Using the paragon research 
engine, we identified >900 proteins per single run. A 
set of 393 proteins met the stringent criteria for being 
quantified, i.e. that a minimum of three specific peptides 
were identities in all 19 shotgun runs. The proteome 
profiles clearly clustered dependent on the PK expression 
level (Fig 1C) This clear separation was caused by 7.8 
percent (31) proteins that were significantly responding to 
changing PK levels. 

These numbers illustrate both the advantages and 
disadvantages of current label free proteomic approaches. 
On the plus side, label-free profiling is much more flexible 
and convenient, as isotope labeled standards are costly, and 
these standards can only be created for a limited number 
of species. In this particular case, i.e. SILAC labeling was 
difficult, as the yeast genetic modifications necessary [44] 
would impact yeast amino acid transport and metabolism, 
and thus potentially interfere with the biological function 
of PK. On the down side, at the current state of the art, the 
number of proteins quantifiable in a labeled experiments 
is higher, as the peptide pools can easily be pre-separated 
before the LC runs. However, also when using isotope 
labels, the number of replicate injections may be limited 
by the throughput of the nanoLC gradients (i.e. the 19 
injections using at a 3hr gradient required one week of 
instrument time). These current bottlenecks for label-free 
approaches might however improve in the near future, as 
recent developments in data-independent acquisition (i.e. 
SWATH [45]) will facilitate significant faster runtime 
while increasing in the number of quantifiable peptides 
and proteins. 

To verify the proteome profiles, we applied the 
targeted MS technique selective reaction monitoring 
(SRM) and analyzed the expression of nine randomly 



chosen PYK targets. Facilitated by the shorter runtimes 
and the high reproducibility of the SRM experiments, 
we included three additional strains that have different 
PK activities, the wild type parent BY4741, as well as 
two yeast strains expressing the PYK2 both at high and 
low level [17]. For seven of the nine tested proteins, 
this extended study did confirm Pyklp, and also Pyk2p 
dependent regulation. Furthermore, the SRM experiments 
allowed grouping of these PK targets according to 
their co-expression pattern. Remarkably, co-expressed 
proteins formed functional clusters: proteins up-regulated 
were glycolytic enzymes, proteins downregulated were 
involved in translation. PK dependence of glycolytic 
enzyme expression is an interesting parallel to the humans 
in respect to cancer development: hydroxylated PKM2 
leads to an upregulation of glycolytic enzymes through 
interacting with hypoxia inducible factor -la (HIFla) 
[38, 39]. Thus, in both species regulation of glycolytic 
enzyme abundances can be attributed to PK, and hence, 
at least also parts of non-catabolic functions of PK appear 
to be conserved between yeast and humans. Interestingly, 
we also observed overlapping targets upon an activation 
of the PPP during the oxidative stress response [16]. As 
low PYK activity increases the concentration of PPP 
metabolites [17], it is conceivable that the PPP is involved 
in gene expression regulation mediated by PK. In this 
context is has to be noted that switches of glycolysis to the 
PPP can howeveronly partly attributed to transcriptional 
regulation, i.e. metabolic shifts in glycolysis and the PPP 
can be regulated on the metabolic level only [46], and may 
occur much faster than transcriptional regulation [47]. 

The vast majority of identified PK targets were 
involved in metabolic processes related to glycolysis, 
the synthesis of nucleic acids, amino acids and proteins. 
Although it is assumed that PK mediates the regulation of 
glycolysis [39, 40], similar mechanisms might also exist 
for amino acid metabolism. To obtain a metabolic picture 
of the global PK function in this process we thus set out to 
use ion exchange chromatography and ninhydrin detection 
to profile the levels of free amino acids in the PK model 
strains. Comparing yeast with different PK expression 
level and activity, we found a clear reconfiguration 
of the overall amino acid profile. In cells with low PK 
activity, seven of 17 amino acids were significantly less 
concentrated, and two were higher concentrated. 

Amino acid levels are important indicators 
for the metabolic homeostasis of cells - and closely 
interconnected with protein metabolism. In yeast, 
mammalian cells, and flies, proteasome inhibition is lethal 
due by causing amino acid shortage [48], implying that 
maintaining sufficient amino acid supply is crucial for 
proliferating cells and tumors. One of the amino acids we 
found at decreased concentration was serine. Remarkably, 
it has recently been described for human cells that PKM2 
stimulates de novo serine biosynthesis [49], and that in 
turn serine can act as allosteric activator of PKM2 [50]. 
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This implies a feedback circle involving PYK that controls 
serine production in mammalian cells, providing evidence 
that the PK's role in regulation of serine metabolism is 
conserved between yeast and mammals. 

Only two amino acids increased in concentration 
when there was low PYK activity, glutamine and 
glutamate. Turnover of these two amino acids may be of 
crucial importance for cancer metabolism, as cells can 
obtain a substantial amount of their energy through the 
process of glutaminolysis [3, 51]. A change in pyruvate 
kinase expression seem thus shift energy metabolism 
towards glutaminolysis, providing sufficient energy supply 
for cancer cells. 

To test for conservation in humans, we searched 
for human orthologues of the proteins regulated by yeast 
PK. Remarkably, two direct orthologues were known to 
be associated with cancer progression. The first protein, 
asparagine synthase, is the orthologue of yeast protein 
ASN1, ASNS (Fig 3). ASNS is as potential target for 
the treatment of castration-resistant stage of prostate 
cancer [52]. Moreover, the expression level of this 
enzyme is an indicator for the efficiency of the enzyme 
drug L-Asparaginase (L-ASP) in ovarian cancer [53]. 
The second protein, serine hydroxymethylstransferase 
(SHMT1), is the human orthologue of the yeast protein 
SHM1. SHMT1 catalyzes the formation of glycine 
and methylene tetrahydrofolate which are important 
precursors for nucleotide biosynthesis, and is discussed as 
a therapeutic target since the late 1980s [54]. 

Taken together, these observations implicated that 
the function of PK in regulating amino acid metabolism 
is conserved between yeast and human. Therefore, we 
continued our investigations by quantifying further 
potential PK targets in a human tumor. For two reasons, 
we selected thyroid follicular adenoma for these 
investigations: First, these benign cancers are characterized 
by a clear upregulation of the PKM2 enzyme [20]. Indeed, 
normal and adenoma biopsies are clearly distinguished by 
the PKM2 level ([20], graphically illustrated in Fig 4A). 
Also PKM1 was expressed in these tumor tissues, but 
at around a 50 fold lower level. In difference to PKM2, 
PKM1 concentrations did not distinguish healthy thyroid 
from follicular adenoma tissue. Second, as adenomas are 
slow growing tumors, one can exclude unspecific effects 
which may result from altered metabolism of the rapidly 
proliferating cells. Finally, the activity of the respiratory 
chain complexes I - IV, determined in both the healthy 
as well as the adenoma biopsies, indicated significantly 
lowered activity of mitochondrial metabolism (Fig 4B). 
Thus, being characterized by increased PKM2 levels and 
a decrease in respiratory chain efficiency, thyroid follicular 
adenoma represented tissue exhibiting classic features of 
the 'Warburg effect' 

SRM assays were tested for all five human proteins, 
but only SHMT1 and the aminopeptidase AAP-S were 
detected at sufficient intensity for reliable relative 



quantification by our LC-MS/MS workflow. Remarkably, 
AAP-S and SHMT1 were both found clearly upregulated 
in these adenoma samples. At the measure of the relative 
expression level, AAP-S levels clearly distinguished 
adenoma and that of control tissue (p < 0.01). Plotted 
against the absolute concentration of PKM2, both targets 
created clusters that distinguished the adenomas from 
the control tissue. Thus, expression of two direct protein 
orthologues to the yeast PK targets correlate with PKM2 
expression in thyroid follicular adenoma. 

CONCLUSIONS 

Pyruvate kinase emerges as a central hub in the 
regulation of metabolic pathways important for cell 
proliferation and cancer. As a glycolytic enzyme, it is 
responsible for the net ATP yield of glycolysis, and the 
control of the redox state. Here, a combined analysis of 
label-free shotgun proteomics, targeted proteomics and 
amino acid profiling reveals that this protein is further 
important for the coordination of biosynthesis of an 
orchestra of amino acids. Glutamine and glutamate 
levels were highly increased in cells with low PK 
expression levels, which seem to connect the process of 
glutaminolysis with the regulation of carbon metabolism. 
Seven of the remaining amino acids were decreased 
in concentration, including serine, which has recently 
been reported to be a direct regulator of PKM2 activity 
in humans [49, 50]. Finally, studies in thyroid follicular 
adenoma indicated that part of this PK function is 
conserved in man, and revealed two protein biomarkers 
of this disease state. AAP-S and SHMT1 correlated with 
the PKM2 expression level and were upregulated in the 
tumor. Thus, PK plays an evolutionary conserved role in 
the coordination of carbon and amino acid homeostasis, 
making it and its downstream effectors promising 
candidates for modulating (cancer) metabolism for 
therapeutic purposes. 

METHODS 



Yeast strains and growth conditions 

ApyklApyk2 yeast strains expressing Pyklp or 
Pyk2p under control of the TEF1 or CYC I promoter were 
described previously [17]. The Apyk2 strain was generated 
from BY4741 by single gene replacement with kanMX4, 
and verified by PCR and SRM. Yeast was grown in 
yeast peptone (YP) media containing 2% dextrose, 2% 
galactose, 3% ethanol/ 0.1% glucose or 2% raffmose as 
indicated. All measurements were performed on mid-log 
batch cultures. 
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Enzyme assays 

Follicular thyroid adenoma and normal adjacent 
thyroid tissues (20-100 mg) were homogenized with 
a tissue disintegrator (Ultraturrax, IKA, Staufen, 
Germany) in extraction buffer (20 mM Tris-HCl, pH 
7.6, 250 mM sucrose, 40 mM KC1, 2 mM EGTA) and 
finally homogenized with a motor-driven Teflon-glass 
homogenizer (Potter S, Braun, Melsungen, Germany). 
The homogenate was centrifuged at 600g for 10 min at 
4°C. The supernatant (600g homogenate) containing 
the mitochondrial fraction was used for measurement of 
enzyme activities. OXPHOS enzyme and citrate synthase 
activity [55-59] were measured as described earlier. 
Activity of respiratory chain complexes is given relative 
to citrate synthase activity. 

Shotgun proteomics 

An TripleTOF™ 5600 system (AB SCIEX, 
Concord, Ontario, Canada) was coupled with a Dionex 
Ultimate 3000 RSLCnano system (Dionex, Sunnyvale, 
California). Details to instrumental parameters and data 
processing can be found in the supplementary material. 
Briefly the raw data was analyzed using ProteinPilot™ 
Software v.4.0 (AB SCIEX) on database UniprotKB/ 
Swiss-Prot 57.15 in thorough search mode applying the 
integrated false discovery rate analysis tool (PSPEP) [60]. 

Subsequent quantification of peptides and proteins 
was performed automatically in Peak View™ Software for 
proteins with <1% false discovery rate as performed using 
PSPEP. Areas of peptides and proteins were automatically 
transferred to MarkerView™ Software for visualization of 
differences for technical replicates and statistical analysis 
using principal component analysis followed by principal 
component variable grouping to identify protein groups 
with similar trends in regulation. 

Selective reaction monitoring 



energy (CE) settings are given in the Suppl. Table 3. 
Amino acid Analysis 

Amino acid analysis was performed using ion- 
exchange chromatography on a Biochrom 30 amino acid 
analyzer with ninhydrin detection. Each sample was 
analyzed in triplicate. 200 (xl sample was added to 200 
(0.1 of 0.5 mM DL-2,4-diamino-n-butyratedihydrochloride 
(used as internal standard) and 4.9 % 5-sulfosalicylic 
acid to precipitate the protein. Samples were mixed and 
centrifuged for 10 minutes at 3345g and the supernatant 
was filtered (Ultrafree-MC Durapore PVDF, 0.22 jim; 
Waters, UK) for 10 minutes at 3345 g to remove all protein 
precipitate. A series of lithium citrate buffer solutions 
were run through a lithium column (8 |am Biochrom, 
UK) containing the samples. Individual amino acids were 
eluted according to their pH. Reaction with ninhydrin 
was utilized to elicit a spectrum of colors at different 
wavelengths giving a chromatogram representing the 
concentrations of each amino acid. 
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